o^-, Scattering in one-dimensional heterostructures 
described by the Dirac equation 



(N 
C 

a 
>—> 

o 



N. M. R. Peres 

Center of Physics and Department of Physics, University of Minho, P-4710-057, 
Braga, Portugal 



Abstract. We consider electronic transport accross one-dimensional heterostructures 
described by the Dirac equation. We discuss the cases where both the velocity and 
the mass are position dependent. We show how to generalize the Dirac Hamiltonian in 
(ZJ | order to obtain a Hermitian problem for spatial dependent velocity. We solve exactly 

the case where the position dependence of both velocity and mass is linear. In the 
case of velocity profiles, it is shown that there is no backscattering of Dirac electrons. 
In the case of the mass profile backscattering exists. In this case, it is shown that the 
linear mass profile induces less backscattering than the abrupt step-like profile. Our 
results are a first step to the study of similar problems in graphene. 

O . 

O , PACS numbers: 72.10.-d, 73.21.Hb, 73.23.Ad 



j> ■ 1. Introduction 

J> 

Most of the accumulated knowledge about the physics of heterostructures assumes that 
the electrons in these materials are effectively described by the Shrodinger equation, 
with a position dependent mass. The common belief is that Dirac equation is of no use 
Q\ \ in condensed matter physics (spin-orbit coupling can be treated using the Pauli version 
of the Schrodinger equation). The discovery of graphene [HE] changed drastically this 
perspective. Condensed matter physicists are now facing a condensed matter system 
where the effective low-energy model for the quasiparticles is that of a ultrarelativistic, 
i. e. massless, Dirac equation, albeit with an effective Fermi velocity that is much 
lower than the velocity of light. In graphene the Fermi velocity is vp = c/300, with c 
the velocity of light. In fact, since the isolation of graphene crystallites [H EI [3] that a 
renovated interest in the properties of the massless Dirac equation in 2+1 dimensions 
started to emerge (3J 0]. Although our interest in this paper is on the backscattering 
properties of electrons effectively described by the Dirac equation in 1+1 dimensions, 
as an effective low-energy theory of the electronic properties of a quasi-one-dimensional 
solid, we shall revise some of the properties of Dirac electrons in the context of graphene, 
given the accumulated amount of evidence for Dirac fermions in graphene. 

How can the low-energy description of a solid be given by an effective Dirac-like 
equation? In one-dimensional physics, this takes place due to the linearization of the 
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spectrum close to the Fermi momentum, which introduces right- and left-movers with 
a linear energy-momentum relation. Another possibility, that we discuss in the bulk of 
the paper, is by having two atoms per unit cell with special type of hopping. In two- 
dimensions, this possibility has been realized by graphene. The electrons in graphene 
are confined to move within the tt— orbitals of system, formed from the overlap of the 
2p z atomic orbitals of a single carbon atom. The transport and optical properties of 
graphene are determined essentially by the behavior of the n— electrons. The electronic 
density is such that graphene has one 7r— electron per carbon atom, and is therefore a 
half filled system. It turns out that the band structure of graphene around the Fermi 
surface is highly non-standard, since the relation between the momentum and the energy 
is linear, with a proportionality coefficient given by Vp ~ 1 x 10 6 m/s and termed Fermi 
velocity. In fact, the Fermi surface is reduced to two points (the so called Dirac points 
K and K') in the Brillouin zone, with a zero density of states at the Fermi-points. This 
fact gives graphene its non-conventional properties. 

Starting from a k ■ p approach [5], adapted to degenerate bands off the T— point 
[6], we obtain an effective Hamiltonian in real space, valid around the K point, having 
the form 

H = v F a-p, (1) 

where cr = (a x , a y ) is a vector using the Pauli matrices and p = —ih(d x ,d y ). The 
Hamiltonian (pp) is nothing but the Dirac Hamiltonian for massless particles in 2+1 
dimensions, with an effective light-velocity vp. Clearly, one expects that the physical 
properties of such a system will be different from those where the Schrodinger equation 
is valid. Note that this situation is an example of complex emergent behavior p2], since 
the original problem was that of independent electrons (and therefore Schrodinger like) 
in a periodic potential. 

Those electrons in graphene lying close to the Fermi points (K and K'), will have 
a different response to external potentials than those in other materials described by 
the Schrodinger equation. In fact, for Schrodinger electrons, if we consider that the 
potential varies at a given point in space from zero to a finite value Vq (the so called 
step potential) there is always a finite fraction of impinging electrons that are reflected 
back. More over, if Vq is greater than the energy of the impinging particles, there will 
be an exponential attenuated wave function in the region where the potential is finite. 

For Dirac electrons the situation is different. It was shown by Klein [16] that Dirac 
particles can not be confined by an arbitrary large potential Vq. In fact, they pass 
through strong repulsive potential without the exponential decay that characterizes 
Schrodinger particles. This is called Klein tunneling. If the particles are massless 
and moving in one-dimension the situation is even more dramatic, since there is no 
backscattered probability flux, no matter how large the potential is fl3l [T4]. 

In addition to external potentials, a particle propagating in solid state systems can 
face a situation where its effective mass changes in space. This possibility takes place in 
heterostructures, where the mass of the particle is written as m(r), r being the position 
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of the particle. In this situation, the Schfodinger equation has to be modified, in other 
to comply with hermiticity and flux conservation. 

If the particles are described by the Dirac equation, the information on the material 
properties is encoded both in the Fermi velocity [6j and in the mass. We can then imagine 
the possibility of producing a herestructure where the Fermi velocity, Vp(r), and the 
mass change with the position of the particle. As we show bellow, even in this case, there 
will be not backscattering if the particle is massless and the velocity is allowed to vary 
from point to point. In the context of graphene physics, a position-dependent Fermi 
velocity was already considered in the case of curved graphene [7]. Another possibility 
of producing both a Fermi velocity and a mass term that are position dependent is to 
subject the system to strain. It has already been shown using ab-inito calculations, in 
the context of graphene physics, that strain leads to gap formation 0, [9]. It is to be 
expected that for general strain, different parts of the material presents local values 
of Fermi velocity and energy gap. Also in the case of graphene, it was shown that 
depositing the material on top of Boron Nitride leads to gap-formation [10J. 

In this paper we discuss several possibilities for the scattering of Dirac particles 
through a region where both the velocity and the mass are position dependent. Although 
our motivation is to study graphene strips [T5j in a latter publication, we start by giving 
two exact solutions for one-dimensional systems. The generalization to a quasi-one- 
dimensional system, such as a narrow nano-wire, is easy to obtained and will be given 
in a follow up publication. In order to see the differences between the Schrodinger and 
the Dirac problems, we start by revising the case of one-dimensional heterosctructures 
described by the Schrodinger equations and latter move the study of the Dirac case. 



2. Warming up: the Schrodinger electrons in ID 

In order to compare how the scattering of Schrodinger electrons differs from that of 
Dirac particles, we give a brief account of the scattering of electrons by the interface of 
a heterostructure. 

Let us consider the case of ID Schrodinger electrons with position dependent 
effective mass m(x), such that m(x) = m~ , for x < and m{x) = m + , for x > 0. 
The boundary conditions in this case are (see final paragraph of this section; for more 
general boundary conditions see Ref. [H] ) 

Vr(o) =^ + (o), (2) 

— ff(0)^^ + (0). (3) 
m ax m + ax 

We consider that the particles are moving from the left to the right. Therefore the wave 

function is 



if) (x) 
iP+( x ) = te ipx 



e ikx _|_ re ikx ^ 



(4) 
(5) 
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where the values of k and p are fixed by energy conservation: 

k2t,2 fc2 2 

E= — = ^. (6) 

2m- 2m+ v ' 

The reflectance \r\ 2 is determined from the boundary conditions as 

,o (m + k — m-p) 2 , . 

\r\ = — , (7) 

(m + k + m~p) 2 ' 

which has the limiting value \r\ 2 — > 1 for m + — > oo. The transmittance coefficient |t| 2 is 

up A(m + k) 2 

[m + k + m~p) 2 ' 

We can check now that the boundary conditions used for the wave function are the 
correct ones, since they satisfy the flux current conservation. Using absolute values, we 
have indeed 

jinc _ jreft _|_ jtrans ^ 

with j mc = hk/m + ', j re ^ = |r| 2 M/m + , and j trans = \t\ 2 Hp/m~. 

The problem just described represents an heterostructure, where two materials are 
glued together (due to similar lattice constants), having different effective masses in the 
two materials. The used boundary conditions are obtained writing the Hamiltonian as 

h 2 d 1 d . . 

H= : z — -— — . 10 

2mdxm(x)dx 

This choice defines an Hermitian problem [j] since J dx(Hip)*ip = J dxijj*Hijj. The 
construction of the ( fTOl) guaranties that the total probability flux is conserved. The 
boundary condition j3j), follows immediately from the writting the Hamiltonian in the 
Sturm-Liouville form. 



3. Dirac electrons in ID 

Let us now study the case of the Dirac Hamiltonian. We are interested in cases where 
both the velocity and the mass of the particles depend on their position in space. We 
study first the case of massless Dirac particles, like those present on graphene. Latter 
we add up the mass term. The Hamiltonian for a massless quasiparticle described by 
an effective Dirac equation, with an effective velocity of light Vp, is given by 

h d 

H = v F a x -— , (11) 

i dx 

% The Eq. (Tl0|l is a particular case of a Sturm-Liouville operator, which has the general form: 

--T- (p( x )t~) + l( x )v = ^w(x)y , 
ax \ ax J 

where A is to be determined from the boundary conditions, and p(x), q(x), and w(x) are given functions. 
Sturm-Liouville problems are Hermitian. 
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where a x is the x Pauli matrix. Let us now consider that we make an heterostructure, 
made of two different materials, both of them described by an effective Dirac equation, 
such as (fTTi l. An example could be a strip of the material where part of it is subjected 
to strain and other part is strain free; this leads to different velocities in the two 
parts of the material. This situation calls for a model where Vp is position dependent: 
Vp = vp(x). Although this situation does not make sense in high-energy physics, it is 
quite conceivable in condensed matter physics, since the value of vp is determined by the 
material under consideration - different materials can have different Fermi velocities. 

The trivial replacement vf — *> vp(x) renders the problem non-Hermitian, as can be 
seen by applying the definition given above, when we discussed the Schrodinger case. 
There is however a way out. It is a simple matter to show that the operator 

H= \Zv F (x)a x - — ^vp{x), (12) 



is Hermitian and reduces to Eq. (fill in the particular case vf(x) = vf. Note that the 
derivative will act on the product ^vf(x)iP(x). We stress that we are not studying the 
scattering of relativistic particles in condensed matter systems, but instead describing 
the scattering of particles represented by an effective low energy model that is formally 
equivalent to the Dirac equation (a situation that takes place in graphene). It would 
be interesting to derive Eq. ( fT2i ) from a microscopic Hamiltonian. In Appendix A we 



give a tight-binding model whose effective low-energy theory is given by Eq. (fT2"l) . The 
problem (fT2l) has spinorial wave function of the form ip* = (tp*,^) and the probability 
flux is computed as S x = VF(x)ip^a x ip, as can be shown using the traditional derivation 
of computing the time change of the probability density [17]. 

The eigenproblem Hip = Eip, with H given by Eq. (fT2"l ). corresponds to two coupled 
first-order differential equations of the form 



/ ; , h d \\/ v{x)ip2{x)} „ , . , . , 

VHxj- LV y 1 n =EMx), (13) 



r —-hd[^/v(x)ip 1 (x)} 

V¥): j- 



E*p 2 (x). (14) 

i dx 

It is straightforward to show that the two first-order differential equations given above 
can be put in Sturm-Liouville form: 

d d 

^ W t Fl(l)1 = ^ l(l) ' (15) 

with yi{x) = y/vF{x)ipi{x), p(x) = vp(x), X = E 2 /h 2 , and w(x) = l/v F (x), and ^2 0*0 
obtained from 

H i — 7— r d 



ip 2 {x) = -^V v F(x)-^-[y/v F (x)ip 1 (x)] . (16) 

If the velocity profile changes continuously, one expects that the continuity of the wave 
function should be valid. If the velocity profile changes abruptly at a given point in 
space, say at x — 0, from v(x) = f_, for x < 0, to v(x) = v + , for x > 0, one has to use 
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Eqs. f fT3l) and (fT4j) to derive the boundary conditions. Defining 2/2 0*0 = V / ^0^)^ 2 ( x ) 
Eqs. (fTBl and (TT41) can be written as 

fi G% 2 (x)] _ 

i ■u(x) ' 

= E y 2 (x) 
i dx v(x) 

Integrating Eqs. (TT7l) and (TT8l) using a symmetric infinitesimal interval around x = 
one obtains the condition 

K(O-) = j/i(0 + ) , (19) 

with i = 1,2. Note that condition (fT9l) implies the discontinuity of the wave function. 
Let us now move to the solution of several particular cases. 

3.1. The step-like velocity profile 

Let us consider the case 

v F (x) = Vp6(-x) +v+6(x) . (20) 
The solution of the Dirac equation reads 




^.(X)= ( , ) f "'■'■ + /■( , ) r (21) 

and 

M*)=t( \ Je*", (22) 

with E = Vpqh and = Vpph. The above wave functions were obtained by solving the 
Dirac equation for x ^ 0, where the velocity is constant. Using these solutions and the 
boundary condition (fT9l ) one obtains 

JvZ{l + r) = v /ttjTi 



(23) 

^(1 - r) = ^ , V ; 



which is satisfied only for r = 0. Had we consider the general case of both a velocity 
profile (1201) and a potential profile of the form 

V(s) = Vo0(aO, (24) 
under the condition £7 > Vq, the boundary conditions would still give r = 0. This result 
is called Klein tunnelling [16]. If we consider the same velocity and potential profiles, 
but now take the energy E < Vq, the boundary conditions still gives the result f l2"3"l) . but 
the wave function of the propagating mode for x > is now different and given by 

i>+(x) =t ( j j e~ ipx . (25) 

In this case the result will also be r = 0. The conclusion is that it is not possible to 
backscatterer massless Dirac electrons with a step-like velocity and potential profiles in 
ID. 



Scattering in one- dimensional hetero structures described by the Dirac equation 



7 



3.2. The linear velocity profile with massless particles 

We have seen that an abrupt change of the velocity at the interface produces no reflected 
particles. Let us now study the case of a smooth change in the velocity from Vp to Vp. 
(The result can be guessed from the outset!) To that end, we choose the profile 

v F (x) = Vp 9(-6 -x) + (v + xA) 9(5 - \x\) + vp' 9(x - 5) , (26) 

where we have defined v and A as 

(27) 



v F - v F 



- 2S ■ < 28 » 
For the cases x < —5 (region I) or x > 5 (region III) the Fermi velocity is constant and 
the solution of the Dirac equation is elementary, as we have seen before. The interesting 
case is therefore the region |x| < 8 (region II). In this case we have to use the Dirac 
equation in the form (Tl2l . Explicitly we have to solve the problem Hip = Eijj with H 
given by 

h d 



H = y/v + xA a x — — Vv + xA . (29) 

i dx 

Writting the differential equations satisfied by the spinors, we obtain for the -01 spinor 
(considering the substitution y = y/v + xAtpi) the equation 

-4(x)0-A, F (-)g = e 2 y, (30) 

with e 2 = E 2 /h 2 . Making the replacement 

v + xA „. , . s 

z= — = — = 9{x), 31 
v 

we obtain the result 

2 d2 V . d y . 2 n fao\ 

z dz^ + z T z +vy = ^ (32) 

with v 2 = e 2 / A 2 . Making the additional replacement u> = hiz, Eq. (l32l) is reduced to 
that of the harmonic oscillator T 

|^ + A = 0. (33) 
The general solution of Eq. (1331 is elementary and from it one obtains y(x) given by 
y(x) = A mn[uLe(x)] + B cos[uL9(x)}, (34) 

with 

L9(x) = In _ . (35) 

§ One should note that the solution of Eq. (1321) can also be obtained by noticing that y = z rn is a 
solution if m 2 = —v 2 , leading to y(z) — Az w + Bz~ w . 
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The component ipi of the spinor is obtained from ipi(x) 
obtained using Eq. (fT6l) and reads 



h , — r -dy(x) 

The boundary conditions are if} 1 {—5) 
written as[jj 

e ikS C-\- 

s + /^/4 

C+A- 

= sm[vL6(±5)) , 
= sm[vL9(±8)\ . 
h\/v^ A v 



if, 



with 



V 

s± 
c± 

A+ 



C-h 







-SL 









M 


_ e iqS 




_A+ 


_ e iqS 



5) and V /7 (^) 
\ 



y(x)/y/v^(x), and ^2 is 
(36) 

if) 1 11 (5), and can be 



/ r \ 

B 
\ t J 



( e- ikS \ 


V 



(37) 



(38) 
(39) 

(40) 



iE v 6{±5) ' 

The fraction of reflected flux is given by |r| 2 and the transmitted flux is 1 
Vp\t\ 2 /vp. The explicit evaluation (which is some what lengthy) of the coefficients gives 



\t\ 



F_ 

+ 1 
F 



(41) 
(42) 



and therefore the electrons are totally transmitted across the heterojunction. Of course 
that this result could have been anticipated from the conclusions of Sec. 13.11 since it is 
always possible to represent a well-behaved function by a sum of infinitesimal rectangles. 
(It is however elegant to have an exact solution to a given problem.) 



3. 3. The linear velocity profile with massive particles 

We can now use what we have just learned to a deal with a more general situation 
where the electronic spectrum changes due to a change of the material. The simplest 
case is that where for x = 5 the mass of the quasi-particles jumps from zero to a finite 
value. The connection between the two materials is represented by the linear profile of 
the velocity discussed in Sec. 13.21 The mass profile is 

m(x) = md{x — 5) , (43) 

that is, the change is velocity takes place outside the region of finite mass. Since the 
particle has a finite mass the spectrum changes to E = ±^Jv\q 2 + 
function in this region changes to 



m 2 v\, and the wave 



ifj lu (x) =t 



E 2 — m 2 v\ 
E — mv 2 



-iqx 



(44) 



Note that the dimensions of A and B are \/L/T, and A ± has dimensions of ^JT/L. 
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The overall modification is the replacement of e iq by \/ E 2 — m 2 v \e iq in the third row, 
and of e iq by (E — mv\)e iq in the fourth row. Working out the calculation of |r| 2 we 
obtain 



r 



|2 



-[l + 2e(-e + v / e 2 -l)] ) (45) 



with e = E/(mv\). As expected, the result does not depend on 5 and on Vp, due to 
the reasons discussed in Sec. I3.2L A plot of Eq. (1451) is given in Fig. [TJ The situation 
is same if the mass profile had been chosen as 

m(x) = mQ(x + 5) , (46) 

this is, the linear profit 



Figure 1. Left: plot of Eq. (f45]) as function of e. Right: plot of Eq. lj45j) as function 
of e 



3.4- The linear mass profile 

We have seen above that the linear velocity profile does not contribute to the 
backscattering of Dirac-fermions, even when it is combined with a region where fermions 
are massive. In fact, any velocity profile produces no backscattering in one-dimension. 
All the backscattering comes from the change in the dispersion, due to the presence of 
the mass term. This fact motivates the question of what is form of the wave function 
and of the coefficient |r| 2 if the mass does not change abruptly but in a smooth way? A 
choice for the change of the mass profile, leading to an exact solution, is that described 
by: 

X 

mix) = mOix - 8) +m-9(x)9(8 - x) . (47) 

o 

We approach the solution of the scattering problem by solving first the Dirac 
equation subjected a mass profile m(x) = m|. The method of solution is inspired 
in that used for the 3+1 Dirac equation [18l [19]. In this case the Fermi velocity is 
constant, and therefore we have to solve Eq. ( fill with the additional term a z mv F ^tp: 

J% dip oc 

vfct x -- h a z mv 2 F -tp = Eip . (48) 

i ax o 
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:<T y <J z . (49) 



h v F-T2 = ~l( E ~ m 2 v%x l J8 l ) + a y Kmv F /5\ijj . (50) 



The solution of this problem proceeds in several steps. The first is to operate on the 
left of Eq. (1481) with the operator 
d 

dx 

After some algebra we obtain the following result 

d 2 ^ 
dx 2 

The two components of the spinor are still coupled in Eq. ([5DJ) . In order to decoupled 
them we use the unitary transformation 

(j)=U^, (51) 

*-(;:)■ 

The above transformation changes o~ y to a y = U^a y U , which, still mixing the spinors in 
the Eq. (l53l) below, does it with the same sign (this is a crucial step). After applying 
the unitary transformation we obtain 

h 2 v 2 S^ = -\{E 2 - m 2 v F x 2 /5 2 ) + a y hmvl/5]<P . (53) 
dx z 

We now introduce two new functions defined by F± = 4>\ ± (f>2, where (pi (with i — 1,2) 
are the components of the spinor (J). In terms of the functions F± the eigenproblem takes 
the form: 

- + [m 2 v 4 F x 2 /S 2 - e 2 ± ]F ± = , (54) 

where e 2 ± = E 2 ±hmv F /5. The general solution of Eq. f l54l) is given in terms of parabolic 
cylinder functions D u (x) (see |Appendix B[ ) 

F±(x) = A ± D v ±(V2(b/afl A x) + B±D v ±(iV2(b/ 'a) l ^x) , (55) 

with 

^= -^ + ^4t = z/ -( 1 T1)/2, (56) 
1 e 2 

4= - %= = -z/f - 1 (57) 

E 2 5 
2Hv F m ' 

a = h 2 v 2 F , (59) 
b =v F m 2 /6 2 . (60) 

If we were looking for the solution of the problem valid to all values of x, the functions 
D u (x) with imaginary argument, since are not real and are not normalizable in the 
infinite volume, had to be excluded. Therefore, in the limit x — > ±oo the normalizable 
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solutions are those where F±(x) represents the one-dimensional harmonic oscillator wave 
functions (see Appendix C ). Thefore the solution in the interval x G [— oo, oo] is 

F ± (x) = 2A ± D ut (V2(b/a) 1 / 4 x) . (61) 

This choice guaranties that the have function is well behaved at x — > ±00 when v x is 
equal to a positiver integer (see Appendix C ). Therefore, in the infinite volume, the 
solution of Eq. f l50l ) takes the form 



= I iA +D ; } - 1 ■ < 62 > 

where the argument of the parabolic cylinder functions has been omitted for simplicity 
of writing. We still have to check whether Eq. (l62l ) is solution of Eq. (l48l ). Introducing 
the solution (l62l in Eq. (T48l) we obtain that the wave function is a solution if 

A_ = ^A + , (63) 




A =v F hJ^i^=v F H/(3. (64) 

What we have discuss so far assumes that Eq. (1481) holds for every x. Our interest, 
however, is on the scattering problem of electrons when Eq. (jl8i) holds for x e]0, 8[. In 
this case, the general solution (1551 ) with both real an complex wave functions must be 
used. The strategy is the same used before. We known at the outset that the solution 
of Eq. ( l48l ) is obtained from the solution (l62l ) by fixing the value of a constant, as in 
Eq. f l63l . The final solution is 



ip(x) = A 



D v {x/P)+vh/ED v .. x {x/P) 
iD v {x/P) - ivA/ED^x/P) 



D_ v (ix//3)+iuA/ED- v .. 1 (ix/p) . 
M -iD_ v {ix/P)-vA/ED„ v _ x {ix/P) 1 ' ' ' 

which is valid in the region x g]0, 8[. As in Sec. 13.21 we wave the incoming wave function 
given by Eq. pTl) . for x < 0, and the transmitted one given by Eq. (jUj) , for x > 5. The 
reflected and transmitted amplitudes r and t, respectively, are obtained from imposing 
the continuity of the wave function at x = and x = 5. In Figure [2] we give some 
numerical examples of our calculation. It is clear from Fig. [2] that the smoothness of 
the change of the mass profile leads to a larger transmittance, for a given energy, than 
when the mass-profile changes abruptly. Also, when the change in the value of the 
mass takes place over a relatively large regions, the transmittance shows the presence 
of resonances, as can be seen in panel (c) of Fig. [2l 
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4. Discussion and conclusions 

We have studied several scattering problems using a modified version of the Dirac 
Hamiltonion which incorporates the possibility of a spatial dependent velocity and mass 
terms. Two exact solutions were given. We showed that in the case of a velocity 
profile it is necessary the modification of the original Dirac Hamiltonian, in order to 
have a Hermitian problem. We have shown that Klein tunneling is not suppressed by 
a change on the velocity profile, with a transmittance equal to unity always. This was 
understood by studying the case of an abrupt change in the velocity profile and also by 
solving exactly the case where the velocity changes linearly across a given region. 

We have also studied the case where the mass term depends on position. For this 
situation, we solved the cases of an abrupt change of the mass value and of a linear 
change of the mass profile. In both cases we see the presence of backscattering, with 
values of the transmittance smaller than unity. The smoother mass profile induces less 
amount of backscattering. 

It is interesting to consider next the case of Dirac electrons in a strip of finite width 
W, a situation relevant for graphene strips, and see if in this case a position dependent 
velocity profile does produce backscattering. This will be considered in a follow up 
publication. 
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Figure 2. Transmittance of Dirac electrons across a region where the mass profile 
changes according to Eq. (|47|) . The dotted line represents the transmittance in the 
case of an abrupt change of the mass profile, as given by Eq. (|45l) . The energy gap is 
represented by the value of mVp. In panel (a) we give the results for a value of 5 of 1 
nm. In panels (b) and (c) we plot the same but for S — 10 nm. In panel (c) we give 
a zoom in of the solid curve of panel (b); it is clear the existence of resonances in the 
transmittance. The Fermi velocity is 10 6 m/s, leading to vfTi = 0.66 eV-nm. 
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Appendix A. A derivation of the Dirac equation with a position dependent 
velocity 

Consider a one-dimensional tight-binding Hamiltonian with two atoms per unit cell such 
the the hopping parameter within a unit cell n is —t n and among nearest unit cells in 
t n . The Hamiltonian in second quantization reads 

H = 22[-t n (alb n + b ] n a n ) + t n {a) n b n - X + 6^_ia n )] . (A.l) 



The in the case t r 
electrons is 



t (the hopping is independent of the unit cell) the spectrum of the 



E a= ± = a2t\ sin{kc/2) | , (A.2) 

with c is the length of the unit cell vector. Close to k = (zero energy) the spectrum is 
linear in momentum and a massless Dirac spectrum is generated with a Fermi velocity 
given by Vf- Let us now assume the general case of a t n dependent on the unit cell 
position and obtain from the Hamiltonian (I A. If) the effective field-theoretical model 
that describes the system at low energy {E ~ 0). Since the momentum close to which 
the Dirac spectrum develops is k = 0, we can write immediately the effective field 
theoretical model as 



H= — 



dyt(y) 



a (y) c — — + c J~. a \V) 



(A.3) 



dy dy 

where we have used the expansion W{y — c) ~ tf(y) + cd[tf(y)]/dy. Integrating by parts 
the second term of Eq. (IA.3I) and using the Pauli matrices to help condensing the results 
we obtain 



H=- c Jdy 

with v(y) = ct(y)/h and W 
Hamiltonian (TT2D. 



(A.4) 



(y)b'(y)]- The first term in Eq. (IA.4I ) is our proposed 



Appendix B. Weber's differential equation 

We give here some basic information on the Weber's differential equation , aiming to 
give the text a selfcontain nature and to fix notation and definitions. Weber's differential 
equation is defined as 

y"{z) + {v + 1/2 - z 2 /A)y(z) = , (B.l) 

and its two independent solutions are the parabolic cylinder functions y(z) = D v (z) and 
y(z) = D_ u _ 1 (iz). Equation f|54l) is of the general form 

y"{x) + {-ax 2 + c)y{x) = . (B.2) 



Scattering in one- dimensional hetero structures described by the Dirac equation 



14 



Making the transformation x = (3z, with (3 given by (3 = (4a) 1//4 and v = —1/2 + c/3 2 , 
we reduce Eq. (IB.2I) to Weber's equation (1B.1I) . Using for a and c the particular values 



of our problem we obtain 



hv F 8 
■) 2 F i 



P = \h^—, ( B - 3 ) 
w 2vim 



and 



1 e\8 , . N 

i/= — + . B.4 

The derivative of the parabolic cylinder functions obeys 

D , v (z)+zD„(z)/2-vD v - l (z) = 0. (B.5) 

D v+1 (z) - zD u {z) + uD v ^(z) = . (B.6) 

Using the results of Ref. [[20]] , the solution of D v (z) can be written in terms of the 
Kummer confluent hypergeometric function, U (a, b, 

D v {z) = 2»/ 2 e- z2 / 4 U{-v/2, 1/2, z 2 /2) . (B.7) 

The Kummer function U(a,l/2,z) is computed using the the Kummer confluent 
hypergeometric function, M(a, b, x), as 

. , ^M(a, 1/2, z) . M(a + 1/2, 3/2, z) , . 

Appendix C. Eigenvalues of the scalar potential V(x) = a z mv 2 F x/8 

The problem we introduced in Sec. 13.41 was that of a particle that moves in a 
heterosctructure with a mass dependent position. We can, however, think of this 
problem as that of a Schrodinger particle moving in the scalar potential V(x) = 
m 2 v F (x/8) 2 . Since Eq. (1541 ) is that of an one-dimensional harmonic oscillator, the 
normalizable solutions have the well known form 

z n (u) = ^Le- u2 / 2 H n {u) , (C.l) 
V 2 n n\ 



where 



Vim ^ 1/2 



« = *{m') (a2) 

upon the identification 

— ~2v 2 F , (C.3) 
m 

(J 2 ^Av F m 2 /5 2 , (C.4) 

with m and co the mass and the frequency of the oscillator, respectively. The spectrum 
is obtained from e 2 ± = 2hv F m{n + 1/2)/ 8, with n — 0,1,2, . . ., leading to 

E\ = 2hv F mn8- 1 , (C.5) 

E 2 _ = 2hv F m(n + l^- 1 . (C.6) 
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Since the energy has to be same for both F + and F_, we choose the solutions = z n 

and F™ = z n _i, with Z-\ = 0. Finally, from the definition F± = (pi ± 02, we obtain 
(properly normalized) 

<j)i{x) = [z n {x) + z n ^(x)]/V2, (C.7) 

fa(x) = [z n (x) - z n ^(x)]/V2. (C.8) 

The above equations are the solution of Eq. (l50l) . 
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If It is interesting to note that the system supports a zero energy mode:F" = z and F_ = 0. 



